Analytical approximations to the core radius and energy of magnetic vortex in thin 

ferromagnetic disks. 



(N 

o 

< 



-a 

o 
o 



>: 

oo 
ov 

(N; 

06 • 

o: 

(N 



X 



Konstantin L. Metlov 
Donetsk Institute of Physics and Technology NAS, Donetsk, Ukraine 8311^ 
(Dated: August 16, 2012) 

The energy of magnetic vortex core and its equilibrium radius in thin circular cylinder were first 
presented by N.A. Usov and S.E. Peschany in 1994. Yet, the magnetostatic function, entering the 
energy expression, is hard to evaluate and approximate. In this communication precise and explicit 
analytical approximations to this function (as well as equilibrium vortex core radius and energy) are 
derived in terms of elementary functions. Also, several simplifying approximations to the magnetic 
Hamiltonian and their impact on theoretical stability of magnetic vortex state are discussed. 
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The first topological soliton was discovered as a so- 
lution of non-linear field theory equations by Anthony 
SkyrmejA It had a form of three-dimensional hedgehog 
and was named subsequently "skyrmion" in honor of the 
discoverer. After the landmark work of A. A. Belavin 
and A.M. Polyakov — topological solitons have crossed 
the boundary into condensed matter physics. The latter 
authors discovered much more topological soliton solu- 
tions in the infinite Hciscnberg ferromagnet, as many as 
there are rational functions of complex variable, map- 
ping any such function into some equilibrium magnetic 
structure. Zeros of numerator (possibly with high mul- 
tiplicity), of these rational functions, correspond to the 
centers of magnetic vortices, while zeros of denomina- 
tor to the centers of magnetic anti- vortices. Thanks to 
the infinite size of the 2D ferromagnet, considered by 
Belavin and Polyakov, solitons appeared as absolutely 
stable. Once magnetic texture with a certain topological 
charge (number of magnetic vortices minus the number 
of magnetic anti- vortices, including the multiplicities of 
corresponding zeros into the count) have been created, all 
the structures with different topological charge appear to 
be separated by an infinite energy barrier. 

The model of Belavin and Polyakov became known 
back in particle physics as non-linear 0(3) a model in 
3+1 dimensions and reformulated elegantly in terms of 
functions of complex variable by G. Woo 4 David J. 
Gross found additional family of "meron" solutions to 
it.— Since then, the original Belavin- Polyakov solutions 
became known as just "solitons" . Merons (and all other 
0(3) a model solutions besides solitons^) have infinite en- 
ergy in unbounded 2-d ferromagnet, but can be realized 
when the ferromagnet is finite.— 

While these solutions were obtained long ago, the ques- 
tion of their stability has a history of its own. Kosterlitz 
and Thouless^ analyzed the stability of vortices in 2D fer- 
romagnet and came to conclusion that they are unstable 
and such order could not exist. It is, indeed, true that 
the energy of Belavin- Polyakov solitons is scale-invariant 
and their size is, thus, undefined. In real ferromagnets, 
however, there are various other interactions (not exotic 
at all), which make the vortices stable. N.A. Usov and 



S.E. Peschany, were the first to show that dipolar mag- 
netostatic interaction fully stabilizes magnetic vortex in 
ferromagnetic cylinder both with respect to core radius 
change 7 - and vortex center displacement £ 

Here, starting from more recent (and more general) 
description of magnetization distributions in finite nano- 
elements via functions of complex variable^ the impact of 
various approximations on vortex stability is reviewed in 
unified manner and the expression for vortex core radius 
in circular cylinder— is re-derived. This expression defines 
the core radius implicitly via an equation and an integral 
of certain special functions, which is very inconvenient 
to evaluate and approximate at small cylinder thickness 
due to analyticity problems. It is, however, possible to 
introduce small parameters and expand the integrals in- 
volved and the vortex core radius into series, deriving an 
explicit analytical approximate (but very precise) expres- 
sions, presented at the end. 

In finite planar nano-elements the equilibrium mag- 
netization configurations can be described via rational 
functions of complex variable with real coefficients^ (as 
opposed to complex coefficients in the case of infinite 
film£i£). The simplest ansatz for magnetic vortex in cir- 
cular cylinder (of thickness Lz and radius R) can be writ- 
ten in this language as 



f(z) = i(z - a)/R v 



(1) 



where z = X + 1Y with X and Y being the Cartesian co- 
ordinates in the cylinder's plane (the magnetization dis- 
tribution is assumed to be independent on out-of-plane 
coordinate Z), Ry is the vortex core radius and a is the 
displacement of the vortex from the origin (a = cor- 
responds to the centered vortex). Let us then define a 
complex function 



w(z, z) 



/(*) < 1 

f(z)/\f(z)\ \f(z)\ > 1 



(2) 



where the line over variable denotes complex conjuga- 
tion. The function w is shown to depend explicitly on 
both z and z because it is, in general, not holomorphic. 
It consists of two parts: soliton (where it is analytic and 
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dw(z,z)/dz = 0) and meron (where ww = 1, joined at a 
line (possibly multiply-connected if there are several vor- 
tices or anti- vortices) |/| = 1. The magnetization compo- 
nents, normalized by material's saturation magnetization 
Ms, are then expressed via stereographic projection as 



2w(z, ' 



1 + w(z,z)w(z,z) 
1 — tv(z, ~z)w(z, ~z) 
1 + w(z, ~z)w(z, ~z) 



(3) 
(4) 



If written via magnetization components, this ansatz is 
exactly equivalent the one by Usov and Peschanyi and 
also belongs to the class of trial functions, considered by 
Kosterlitz and Thouless.- Following the latter work, let 
us first take into account only the exchange interaction. 
In complex notation the exchange energy density (omit- 
ting the factor C/2, where C is the exchange stiffness) 
can be directly expressed via the function w: 



i—x,y,z 



dw dw dw dw 



(1 + ww) 2 \ dz dz dz dz 



(5) 



where d/dz = (d/dX - id/dY)/2 and d/dz = (d/dX + 
id/dY)/2. The total exchange energy can be obtained 
by integrating this density over nano- element's volume. 
Recalling the Riemann- Greene theorem 



-/ «(C,C)dC = 

2 « Jdv 



v 



du(z, z) 
dz 



dX dY, 



(6) 



where u is a complex function of the complex argument 
(not necessary analytic^), it is possible to reduce the 
area integral over cylinder's face T> for the total exchange 
energy to a contour integral over its boundary dT>, pro- 
vided there is a complex function, whose derivative over 
z yields the exchange energy density ([5]). Luckily, such 
function (actually two functions, one for soliton and one 
for meron part of w) can be easily obtained by direct 
integration: 



u s (z, z) = — 
u M {z,z) 



1 df 



l + f(z)f(z)f(z) 8z' 
1 df 



f[z) dz 



log {mm) 



(7) 



(8) 



Thus, from ([5]), the total exchange energy inside the soli- 
ton is 



Eix _2 



1 df(Q 
CL z /2 * .7I /(c) |=i /(C) d( 



dC, 



(9) 



where the fact that |/(C)| = 1 on the integration con- 
tour was used and the additional minus sign appears be- 
cause the original contour of integration had to be walked 
clockwise. The function under the integral is analytic ev- 
erywhere except the vortex centers Zi, where f(Zi) = 0. 
Assuming that line |/(C)| = 1 does not cross the particle 
boundary, it is possible to tighten the contours around 



each topological singularity (vortex or anti- vortex center) 
and use the residue theorem 



dm 



CL z /2 



m dz 



(10) 



In particular, for f(z) from Eq. [T] this gives 
E° /(CLz/2) = 47r. For an arbitrary configuration of 
vortices and anti- vortices this number will be multiplied 
by their total number, including multiplicities. 

For the meron part, the integration boundary is 
multiply-connected. However, on the inner boundaries 
(encircling solitons) |/(C)| = 1 and u M (z,~z) ~ logl = 0. 
Thus, only the integral over the cylinder's outer bound- 
ary remains 



E^k __ 1 I 1 df(0 



CL z /2 2»7m,/(C) d( 



log/(C)/(C)dC (11) 



for f(z) from Eq. [T] and the nano-element, shaped as 
circular cylinder (dT> is \z\ = R) 



2tt /i / \ / t->\ i ( a 2 — 2aR cosfo^+i?.' 1 
£JM f (1-0C0S(^)/Je)l0g( — 



dtp 



CLz/2 J 2{a 2 /R 2 - 2acos(ip)/R+ 1) 

o 

^log(l-^)-2,l„g(^), (12, 

and the total exchange energy eEx = 
E M )/(fi "fBMgTTR 2 Lz) (in subsequent text all the di- 
mensionless energies, denoted by small letter e with dif- 
ferent sub-/superscripts, are in this normalization) is 



T 2 
R 2 



'^=^(2-log^ + l«. M /:i-^|. (I:!! 



R 2 



where — 47r, fj,Q = 1 in CGS units and 75 = 1 in 
Slii and the exchange length L E — ^C/(^ B Ml)^. It 
can be seen immediately that this energy decreases with 
increasing the vortex core size Ry . The above expression 
is, formally, valid only for Ry < R, but it can be easily 
shown that the energy continues to decrease for larger 
Rv, reaching equilibrium for Ry — > 00. This confirms 
the conclusion of Kosterlitz and Thouless 6 - that magnetic 
vortices are unstable when only the exchange interaction 
is taken into account. 

There is, however, another interaction, present in all 
magnets. Namely, the long-range dipolar interaction be- 
tween the local magnetic moments. Strictly speaking, 
this interaction is not instantaneous and its speed is lim- 
ited by the speed of light. The account for retardation ef- 
fects, however, contributes to the dissipation^ It is con- 
venient (if magnetic nano-elements are small enough and 
characteristic timescales are large enough) to consider 
the dipolar interaction in magnetostatic approximation, 
making it non-local. This non-locality still poses a ma- 
jor mathematical difficulty, since it makes the equations 
for equilibrium magnetization distribution not only non- 
linear partial differential, but also integral. To alleviate 
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this difficulty a number of local magnetostatic approx- 
imations had been developed. The most common (and 
very useful for considering domain walls in thin films) is 
based on using the local uniaxial in-plane anisotropy term 
Kms^z instead of magnetostatic interaction. Selecting 
the local bulk anisotropy Kus = Mo7bMJ/2, one gets 
the exact correspondence between the approximate and 
exact magnetostatic energy density in two limiting cases: 
when the film is magnetized in-plane (in this case magne- 
tostatic energy is 0) and out-of-plane (in which case the 
energy is ^ 7s^|/2). Nano-elements have additional 
side surfaces and it was recently proposed by Kohn and 
Slastikov 14 to use a similar expression for local surface 
anisotropy of magnetostatic origin on all surfaces with 
some a priori unknown constant K s , replacing the mz 
by a normal magnetization component on the surface. 
Let us try this approach. 

When vortex is completely inside the particle and is 
centered (a = 0) the meron does not contribute to the 
anisotropy energy, and the contribution of soliton part is 



e\" = — 
R? 



(r 2 +R 2 v Y 



■ dr = 



R 2 V (3 -4 log 2) 
R? 2 



(14) 



The total energy density eEX + at a = now has 
a minimum when Ry = Le/\/3 — 4 log 2. Or, approx- 
imately, i?0 « 2.09698LB = 0.59155y/^L E , indepen- 
dent on cylinder's thickness Lz- This thickness inde- 
pendence is the result of expressing the magnetostatic 
energy in the form of surface anisotropy, and, as will 
be seen later, is wrong. Nevertheless, unlike the purely 
exchange approximation, the vortex core size is now sta- 
ble. It is also worth noting that this approximation is 
exact in the limit of vanishing film thickness, so that 
Ry/Le = lim.\^o Pv(A), where pvW is the vortex core 
radius, computed with full treatment of magnetostatics 



But stable vortex core size is not all, the vortex 
must also be stable with respect to the displacement 
of its center. To consider this, accounting of the 
magnetostatically-defined anisotropy on the cylinder's 
face is not enough, since (for the case of vortex inside the 
particle) its energy is independent on the vortex center 
displacement a. But the exchange energy (|13l) decreases 
when vortex is displaced (|a| increases from 0) and this 
leads to instability. To consider this case properly within 
the magnetostatic anisotropy approximation let us follow 
the proposal of Kohn and Slastikov^ and introduce ad- 
ditional surface anisotropy K$ on the cylinder's side. It 
gives the following contribution to the energy of displaced 
vortex, assuming it is fully within the particle 



E 



a = KsL 2 



2tt 

*/(« 



e-^f(Re^) 
\f(R^)\ 



dip 



R 



(15) 

The exchange energy (f]"2"|) can be expanded as -Eex(&) ~ 
const— C L zTra 2 / (2R 2 ) . Equating two second order terms 



in a gives the condition for vortex stability with respect 
to displacement: R > Rg = C/(2Ks). For radii, smaller 
than Rg : the vortex is unstable. This is, again, only 
partially correct. The stability condition turns out to be 
independent on Lz, which means that, while the par- 
ticles of very small radii are correctly single-domain, in 
particles of disappearing thickness the vortex state is un- 
conditionally stable, which is qualitatively wrong. Nev- 
ertheless, if one deals with particles of specific size and 
considers K and Kg as free parameters, the approach of 
Kohn and Slastikov^ may yield a reasonable approxima- 
tion to the stability and evolution of vortex state, in this 
case Kg will have to vanish as Lz — > 0. The advantage 
of this approach is simplicity, as it allows to get explicit 
expressions for most interesting quantities. The full ac- 
count for long-range magnetostatic interaction, which is 
necessary to build the vortex state theory without free 
parameters, is much harder to do. Yet, in the follow- 
ing text, approximate expressions for vortex radius and 
energy with full account of magnetostatics are derived, 
which are almost as simple. 

To compute the magnetostatic energy let us use the 
magnetic charges formalism, introducing a magnetic 
charge density — div fh which is automatically equal to 
the normal component of magnetization on the surface 
of magnetic material (in which case it is a surface charge 
density a). In centered vortex a = there is only a face 
charge (surface charge on cylinder's face, proportional to 
m 2 ), equal to 



(16) 



v 



where r < Rv because all the charge is concentrated in 
the vortex core. The interaction energy of two such sys- 
tems of charge at parallel planes (faces of the cylinder), 
separated by distance h, can be directly written as 

Rv2ttR v 2tt 

4irU(h) f f f f o-(ri)cr(r 2 )ri dri dyir 2 dr 2 d</3 2 



PolB J J J J \Zr 2 +r^-2rir 2 cos(<z>i-<^ 2 ) + /i 2 
oooo v 

(17) 

It is possible to obtain two equivalent representations for 
this integral, one by directly integrating over the angles 



\Z r i+ r 2- c 2t\T2 cos(<^i - ip 2 ) + h 2 



y/h 2 + (r*i + r 2 ) 2 



(18) 



where K{k) is a complete elliptic integral of the first kind, 
which gives 

*{l (1 + Pl)(l + PIUX 2 + {P1+P2) 2 ' 

(19) 
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where the dimensionless quantities u = 
U/^q'jbMI-kR^), x = h/R v , px = rx/Rv, p 2 = r 2 /R v 
have been introduced. Another representation can be 
obtained using summation theorem for Bessel's functions 
of the first kmcUS 



OO 

!■ 




J (kp)pdp\ dk, (20) 



which is shorter on paper, but, unlike (fH)1) . is, actually, 
a triple integral. The magnetostatic energy of the vortex 
core per unit of cylinder's volume is then 



Ems 



R 3 v 



P^IbM^LzR? L z R 2 



e MS = — ; 7 t — 7v> = t n9 ( M (0) _ u(L z /Rv)), 

(21) 

where the first term accounts for the face charge's self 
energy, while the second for interaction of charges on the 
opposite faces. The total dimensionless energy of the 
cylinder with centered vortex is 



7-2 

R 2 



2 -log 



Rv 



eMS- 



(22) 



Minimization of this energy with respect to Ry results 
in the following equation 



1 

Pv 



3p 2 v (u(0)-u(X/p v )) 
X 



+ p v u'(X/p v ) =0, (23) 



where pv = Rv/Le, A = Lz/Le and prime means 
derivative. This equation in independent on particle ra- 
dius. 

It is, of course, possible to evaluate the integrals (fT^jl . 
(f20| on computer (but even this is tricky, since the sec- 
ond is a badly converging oscillating improper integral 
and the first contains a peak at pi — p 2 , turning into a 
line of integrable logarithmic singularities when h = 0) 
and solve the transcendental equation (|23[) numerically, 
but this is far less convenient (and useful), compared to 
having their simple analytical expressions. Let us now 
obtain such expressions approximately. 

The simplest is the case of large cylinder thickness, 
corresponding to x > 1- In this case the outer integral 
in (|20[) is converging very fast, and also the integrand in 
(fl9| is well behaved. This allows to perform straight- 
forward Taylor's expansion of the integrand and perform 
the integration term by term, which gives: 

(log(4)-l) 2 , 3 + 81og 2 (2)-101og(2) 



35 



4 X 

41og(2)(181og(2) 



8x 3 
25) 



32x 5 



(24) 



Solving (I23p with this u results in the following expansion 
for the equilibrium vortex core radius 



A 
3u 



1/3 



(log(4) - lK 7 , (log(4) - 1) 



(25) 



where uq = u(0) = 0.0826762. Using this expansion, 
the equilibrium energy of thick cylinder (A » 1) with 
magnetic vortex can be written as 



oEQ 



7- 



3p 2 



(log(4) - l) 2 
12p 2 (3ugA 2 ) 1 /3' 



(26) 



where p = R/Le- 

These expressions are simple and for A > 2 are precise 
to a few percent (and for X/^/in > 1 the precision is bet- 
ter than 1%). The problem, however, is that assumption 
of uniformity of magnetic texture along Z axis is not a 
good approximation for thick cylinders (A ^ 1), which, 
eventually, start to develop a 3D structure (such as vari- 
ation of vortex radius with Z at first). In other words, 
this approximation is precise mostly in the region, where 
the vortex ([lj is far from the ground state of the system 
(Eq. [26] may still be useful to find the extent of this re- 
gion by comparing it to the energy of other magnetization 
textures). 

It might be tempting to expand the magnetostatic 
function u(x) around x = an d build the approximate 
vortex state theory on top of that. The difficulty is that 
u(x) is not analytic at this point (it has terms, propor- 
tional to x l°g X) an d the integrals for higher order Taylor 
expansion terms do not converge. Since the very thin par- 
ticles are single-domain and also thin cylinders with large 
radius start to develop a domain structure or several vor- 
tices, bound as finite fragments of cross-tie domain walls, 
such approximation would also be the most precise in the 
region, where the vortex is not the ground state. 

This suggests the idea to build the magnetostatic func- 
tion expansion around an intermediate point ^ = 1, 
where the function u(x) is analytic and otherwise well- 
behaved. Such expansion is most precise for Ry ~ Lz, 
where all the physical assumptions of the vortex state 
theory are valid. This region is also close to the triple 
point on the magnetic phase diagram^ The point x — 1 
corresponds to 



An 



1 



vft(u(0)-t*(l))+u'(l) 



2.7284 « 0.76967V47T 



(27) 

The expansions for equilibrium vortex radius and energy 
about the point Ao are the following 



i=0 

EQ _ log(p/Ap ) , 1 



(28) 



, -J>(A-A r, (29) 



where the first few coefficients are 



i 
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0.189400 


-0.012521 


0.001182 


-0.000093 


bi 


2.387556 


-0.082425 


0.010332 


-0.00171 


0.000315 



5 




cylinder thickness, X=L.JL L 



FIG. 1. Vortex core radius and energy density, as function of 
the cylinder thickness. The dotted and dashed lines show ap- 
proximate analytical expressions, the solid lines are the result 
of the exact numerical calculation. 



These coefficients are just universal dimensionless con- 
stants, they decay rapidly and are sufficient to obtain the 
values, precise up to a 0.5% in the interval < A < 6. 
Comparison of the above analytical approximations with 
exact numerical values of core radius and energy are 
shown in Figure. 



CONCLUSIONS 

Starting with the description of magnetization distri- 
butions via complex variable, various simplifying physical 
approximations for magnetic Hamiltonian were consis- 
tently (in the same notations and units) presented and 
compared with their advantages and deficiencies high- 
lighted. The formulas for the exchange energy of such 
distributions (|10[) - (fTTj) , provided the vortices and anti- 
vortices are fully contained inside the particle (which cov- 
ers many distributions of Ref. 0) , were presented here for 
the first time. Two simple and explicit analytical approx- 
imations for equilibrium vortex core radius and energy in 
circular cylinder were derived, which, together, cover the 
whole range of cylinder geometries. 
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